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Abstract 

In the case of a large class of static, spherically symmetric black hole solutions in higher 
order modified gravity models, an expression for the associated energy is proposed and identi- 
fied with a quantity proportional to the constant of integration, which appears in the explicit 
solution. The identification is achieved making use of derivation of the First Law of black 
hole thermodynamics from the equations of motion, evaluating independently the entropy 
via Wald method and the Hawking temperature via quantum mechanical methods in curved 
space-times. Several non trivial examples are discussed, including a new topological higher 
derivative black hole solution, and the proposal is shown to work in all examples considered. 

PACS: 04.50.Kd; 04.70.Dy; 97.60.Lf; 95.30.Sf 

1 Introduction 

Recent observational data imply an accelerating expansion of the visible universe, which gives rise 
to the so called Dark Energy issue. 

There exist several descriptions of this acceleration. Among them, the simplest one consists in 
the introduction of a small positive cosmological constant in the framework of General Relativity 
(GR), the so called A-CDM model. A generalization of this simple modification of GR consists in 
considering modified gravitational theories, in which the action is described by a function F{R) 
of the Ricci scalar R (see for example [HH]). Typically these modified models admit the de Sitter 
space as a solution and the stability of this solution has been investigated in several places (see for 
example [31[11[S1[5]). Furthermore, viable F{R) models, that is the ones which are able to pass the 
local gravitational GR tests, as well as to describe the inflation with dark energy in a unified way, 
have been recently discussed [71 [51 HI [TU] . Another very interesting class of modified gravitational 
models in which the square root of the quadratic Weyl scalar appears have been investigated in 
Ref. HB. 

Static, spherically symmetric solutions have been investigated in several papers (the simplest 
one being the Schwarzschild-de Sitter solution), and they have been discussed for example in 
Refs. |12[ [T51 HM [TT] . Within this class of higher order gravitational models, the issue associated 
with the energy (mass) of black hole solutions is problematic, and several attempts in order to 
find a satisfactory answer to that problem have been proposed (see for example [121 Uni |T71 [THj 
and references therein). 
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To start with, let us remind the case of GR, in which several notions of quasi-local energies 
may be introduced. In particular we mention the so called Misner-Sharp mass, which has the 
important property to be defined for dynamical, spherically symmetric space-time |19) . where the 
use of invariant quantities plays a crucial role j^HlllI]- For the sake of completeness, we recall that 
in four dimensions, any spherically symmetric metric can locally be expressed in the form 

ds^ = -iij{x')dx'dx' -f n^{x')dnl , i, j e {0, 1} , (1) 

where dVl.\ here is the usual metric on the two sphere , but it could be be the metric of a generic 
two-dimensional maximally symmetric space. Of course, in such cases the black hole will have a 
different topology. The two-dimensional metric 

d-f^ =j^j{x')dx'dx^ (2) 

is referred to as the normal one. The related coordinates are {a;*}, while 7l{x^) is the areal radius, 
considered as a scalar field in the two dimensional normal space. A relevant scalar quantity in the 
reduced normal space is 

x{x) = f^{x)dm^)d3n^)^ (3) 

since the dynamical trapping horizon, if it exists, is located in correspondence of 

X{x) = , (4) 

H 

provided that dix\H 7^ 0. (We use the suffix \h for all quantities evaluated on the horizon). The 
quasi-local Misner-Sharp gravitational energy is defined by 

EMs{x)^\n{x)[l-x{x)] ■ (5) 

This is an invariant quantity on the normal space. Note also that, on the horizon, Ems\h = 
^TIh ^ E, E being the energy of black hole. Furthermore, one can introduce the Hayward surface 
gravity associated with this dynamical horizon, which is given by the normal-space scalar 

(6) 



KH = 7;OjT^ 



H 



being the Laplacian corresponding to the 7 metric. In the spherical symmetric, dynamical 
case, it is also possible to introduce the Kodama vector field /C. Given the metric ([J) it is defined 

by 

rix) = e'^djU , /C'' = = /C^ , (7) 

e^^ being the completely antisymmetric Levi-Civita tensor on the normal space. 

Assuming Einstein equations, in a generic four-dimensional spherically symmetric space-time, 
a geometric dynamical identity holds true in general. This can be derived as follows. Let us 
introduce the normal space invariant 

r(2) = 7*^T„ , (8) 

which is the reduced trace of the stress energy tensor T^i,. Then, making use of Einstein equations, 
it is possible to show that, on the dynamical horizon (see for example|tl9j) 

^H = ^ + 2mZHTl^^ . (9) 
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Introducing the horizon area Ah and the (formal) three-volume Vh enclosed by the horizon, with 
their respective "thermodynamical" differentials dAn = STTTZHdlZn , and dYn ~ AnTZ^dTZH (we 
are assuming a horizon with the topology of a sphere) , we get 



!lILdAH = df'^ 



Stt 



V 2 



t: 



(2) 



" dVH. 



(10) 



This equation can be recast in the form of a geometrical identity, once the Misner- Sharp energy 
at the horizon has been introduced. It reads 



(11) 



In the following, we shall restrict the discussion to the static case in the absence of matter. 
This means that we shall consider only vacuum static solutions. In such a case the metric in ([l} 
can be written in the simpler form 



dr^ 



ds" - -B{r)e'^^^Ue + ^ + r^d^ll 
^ ' B{r) ^ 



(12) 



where a{r) and B[r) are functions of r. Of course the general formalism is also valid in the static 
case, and leads to the horizon condition 



The Kodama vector reduces to 



-a{r) 



(13) 
(14) 



When a{r) = 0, which corresponds to case of GR in vacuum, the static Kodama vector coincides 
with the usual Killing vector (1,0), and Hawking temperature of the related black hole reads 



Tk = 



1 dBirn) 
47r dr 



(15) 



This is a well known result, and it can be justified in several ways, for example making use of 
standard derivations of Hawking radiation |22|). or by eliminating the conical singularity in the 
corresponding Euclidean metric, or making use of the tunneling method, recently introduced in 
Refs. [231 [2311^, and discussed in details in several papers. 

However, as we shall see in explicit examples, within modified gravity it happens to deal with 
black hole solutions with a{r) ^ 0. In this case, the Kodama vector does not coincides with the 
Killing vector. Then one may introduce two Hawking temperatures, the Killing temperature (see. 
Appendix I) 



1 



d{e^o,ir)B{r)) dBir) 



Att V dr dr 

and, making use of ^ the Hayward temperature 



1 



47r 



dr 



Th = 



KH _ I dB{r) 



27r 47r dr 



(16) 



(17) 



which is trivially related to the previous one by Tk = e"''""''T^f. If a{r) = we recover Eg. p^ . 
namely Tk — Th- A detailed discussion about this issue can be found in Refs. [ 2111 121| . in which 
also the dynamical case is discussed. 

In the static case, all derivations of Hawking radiation (for example, the tunneling method in 
Appendix I) leads to a semi-classical expression for the black hole radiation rate 



r = e 



(18) 
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in terms of the change Ai?/^ of the KiUing energy Ek [H], but if one uses the Kodama energy 
Eh for the emitted particle, one has 



AS, 



r EE e . (19) 

From the Eqs. (fT8| and (fT9|) one arrives at the identity 

AEh AEk 



Th Tk ' 



(20) 



which may interpreted as the First Law of black hole thermodynamics as soon as F = e""^"^, with 
S the entropy of the black hole itself. As a result, in the static case the two temperatures Tk and 
Th are equivalent. 

With regard to entropy of the black hole, it is well known that in GR the so called Area Law 
is satisfied, and we have 

5w = ^. (21) 

In GR and in the static case, the First Law of black hole thermodynamics in vacuum reduces to 

dE = TndSw , (22) 

where E is the Misner-Sharp energy evaluated on the horizon. 

Now we come to the key point of our proposal. For a generic modified gravity theories, for 
example the F{R) models, where R is the Ricci curvature, it seams very difficult to define in a 
reasonable way the analogue of the local Misner-Sharp mass (see Ref. [H]). As we will see, an 
exception is the higher-dimensional Lovelock gravity |26| . 

For this reason, in this paper, an attempt is made for obtaining an expression of energy 
associated with black holes solutions in higher order modified gravitational models. The proposal 
consists in the identification of the black hole energy with a quantity proportional to the constant 
of integration, which appears in the explicit solution. The identification is achieved making use of 
derivation of the First Law of black hole thermodynamics from the equations of motion, evaluating 
in an independent way the related black hole entropy via Wald method [37] (see the Appendix II) 
and the Hawking temperature via the quantum mechanics in curved space-time, for example the 
tunneling method [23 or other standard equivalent methods. 

This approach is also supported by the results obtained in Refs. |281I29| . where, on quite general 
grounds, generalizing the Jacoboson results on GR (see the seminal paper [30]) i the equations 
of a modified gravitational theories are shown to be equivalent to the First Law of black hole 
thermodynamics. As it is well known, this issue may be of high relevance in substantiating the 
idea that gravitation might be a manifestation of thermodynamics of quantum vacuum |31| . 

The paper is organized as follows. In Section 2, the Lovelock gravity [26J is revisited, and 
the approach here proposed is shown to work for such a case. In Section 3, the four dimensional 
modified gravity models of the F{R) type are investigated, and the method proposed is applied to 
several cases in Sections 4 and 5. In Section 6 a new topological black hole solution is discussed 
and the method is shown to work, as in Section 7, where the conformal Weyl gravity black holes 
are considered. Finally Section 8 contains the conclusions. In two Appendices, for the sake of 
completeness, the tunneling method and Wald entropy method are briefly discussed. 



2 Lovelock Black Hole Solutions 

In this section, as warm up, we review Lovelock theory with the related static and spherically 
symmetric black hole solutions. This theory is a very interesting higher dimensional generalization 
of Einstein gravity. In general, by making use of higher order geometrical invariants in the action, in 
the metric formalism for the field equations one obtains fourth order partial differential equations. 
However, as Lovelock had shown, one can obtain second order differential equation by making use 
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of higher dimensional extended Euler densities, the so called m-th order Lovelock terms defined 
by 



r _ ^ jrAiCTi---A„iCTm p Piki 



m = 1, 2, 3, 



(23) 



where R\a^'^ is the Riemann tensor in arbitrary D-dimensions and Sp^^^y/.p^^^ is the generalized 
totally antisymmetric Kronecker delta defined by 



The action for Lovelock gravitational theory reads 



5'.: ) 



1= d^Xy/^ 



(24) 



where we defined the maximum order k = [{D — l)/2] and are arbitrary constants. Here [z] 
represents the maximum integer satisfying [z] < z. Hereafter we set a\ = 1. 

For such a kind of theory, the equations of motion in vacuum are second order quasi-linear 
partial diff'erential equations in the metric tensor and read 



5/ = 0, 



(25) 



the Lovelock tensor Qp^ being given by 



1 a. 



2m+l 



m 



PlKl 



(26) 



As we said in previous Section, we shall focus our attention on static, spherically symmetric 
solutions, thus we look for metric of the form 



ds^ = -B{r)df 



dr^ 
B{r) 



'dni 



(27) 



where dfi^ is the metric of a n-dimensional sphere S"" {n — D ~ 2), Such kind of theories become 
quite interesting for D > 4, the four-dimensional case being equivalent to Schwarzschild-de Sitter, 
since Ci = R and C2 is equal to the Gauss-Bonnet quadratic term, which in four-dimensions is a 
topological invariant. 

A direct evaluation of field equations gives [32J 



n d[r"+iTy(r)] 



V - 



2r" dr 
1 [r"+iiy(r) 



n-l 



where W is given by 

m=2 

with a„i = ttm Y[l=i^{n - p)- 



2r 



m^—2rn 
I 



d^r 



[1 - B{r)]r- 



2A 



n(n -f 1) 



(28) 
(29) 

(30) 
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For example, for 13 = 4, fc = 1, and so one has the Schwarzschild-de Sitter solution, while for 
D = 5, k — 2, there is one Lovelock non trivial term (the Gauss-Bonnet, which in five-dimensions 
is not a topological invariant) and one has the Boulware-Deser solution [33J. For higher dimensions 
one has an algebraic equation of increasing complexity, but, as we shall see in the following, for 
our purposes it will be not necessary to know explicitly the expression for the solution B(r). 

For the static metric in (P7)l one has the Killing vector = (1,0) and since 

V,^;: = 0, Q^u = Qu^J., (31) 

the vector = Q^i,V'^ is covariantly conserved and gives rise to a Killing conserved charge. This 
corresponds to the quasi-local generalized Misner-Sharp mass which reads 



where E is a spatial volume at fixed time, rfS^ = ((iS,0), and assuming spherical horizons, 

V \iLn) — r(ri/2+l/2) ' 

In the absence of matter Eq. ((28)) can be integrated and one has 

r"+iT4^(r) =C, (33) 

C being a constant of integration which we will show to be related to the mass of the black hole. 
On shell, that is at the horizon r = rn, B{rH) = 0, Eqs. ([5^ and ([55]) leads to 

IdttG 

Now let us show that a First Law of black hole thermodynamics holds true, with the "energy" 
of the black hole solution, namely the Killing charge obtained below, proportional to constant 
of integration C. In the case of Lovelock gravity the validity of the First Law of black hole 
thermodynamics has been investigated in many places (see for example [Ml 1551 [Ml I37| ). For the 
static case we present a direct and simple proof. 

First of all we introduce the horizon defined by the existence of the largest positive root rn of 

Bir,)^0, ^^0. (35) 
Then from Eq. ([55)) we have the identity 

" ^ m " " n(n + l) ^ ' 

m=2 ^ ' 

On the other hand, taking the derivative with respect to r of Eq. (I33p and putting r — rn, and 
making use again of Eq. p3p , we obtain 



^ n- 
^ a„(n + 1 - 2m) ^^i^2m , /„ _ 1 V"-l - ^^^g _ [ \^ ^ ^n+2-2m 



n+1 



771 



m— 2 \m— 2 




Now, let us compute the "thermodynamical" change of C with respect to a small change of th- 
>From Eq.((53]) one has 



\ ^-^ m n 

\m=2 / 

Making use of Eq. ()37p this expression may be rewritten in the form 

= ^ f E ^rar^-'-'"' + dTH . (39) 
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Let us interpret the r.h.s of the latter identity. Here we are deahng with a static, spherically 
symmetric metric admitting a Killing vector. If there is an event horizon located at ru, then the 
Hawking temperature of the related black hole is given by Eq. (|15p . 

Now, all therniodynamical quantities associated with these black holes solutions can be com- 
puted by standard methods. In particular, the entropy can be calculated by the Wald method 
[m [351 [Ml or other methods if you hke, and one has (see for example |55l[55lHD] l 

27ry(r?„) n ( am 2-2rn\ ( .(.^ 

Sw = Q n 1 + " — TT, T^-^H ■ (40 

SttG \ ^-^ n + 2 — 2m j 

\ ra—2 / 

As a result, from Eqs. (|34p . and (I40p . one has the First Law of black hole thermodynamics 
for Lovelock gravity, that is 

Tk dSw = dEK ■ (41) 

We have shown that for a generic Lovelock gravity, the First Law of black hole thermodynamics 
holds and one can identify the energy of a static, spherically symmetric black hole with the constant 
of integration and Killing conserved charge. 

The generalization to topological Lovelock black holes has been investigated in PT] , and again 
the First Law of black hole thermodynamics has been shown to hold. 



3 F{R) four- dimensional modified gravity 

In this Section we will come back to = 4. To begin with, we recall that the action of modified 
F(i?)-theories reads 

1 



d^xy^F{R) , 



(42) 



where g is the determinant of metric tensor g^^, and F{R) is a generic function of the Ricci scalar 
R. For dimensional reason, F{R) may contain a multiplicative functional dependence on G, the 
Newton constant. 

The equations of motion in vacuum for a general F{R) model read 



(43) 



Here, is the Ricci tensor and the part of 'modified gravity' {MG) is formally included into 
the tensor Gjf^p, which is given by 



p, 



(44) 



The prime denotes derivative with respect to the curvature i?, is the covariant derivative 
operator associated with g^^, and □</> = g^^V i,4> is the D'Alembertian of a scalar field 4>. The 
trace of Eq. (|43|l gives 

iUF'{R) + RF'{R) - 2F{R) = , (45) 

which shows that there exists an additive scalar dynamical degree of freedom represented by F'{R). 
In the metric (IT2]) the scalar curvature reads 



R 



B{r) 



dr 



a{r)~2B{r) 



dr 



a (r) 



f^B{r)~2B{r)^a{r) 



-2- 



(46) 
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In Ref.|TF] the following equations of motion have been found by Lagrangian methods ^\ H51 15]: 



e^CO i^RF'{R) - F{R) - 2F'{R) 
+2B{r)F"{R) 



(l- B{r) -r{dB{r)/dr)) 



2 dB{r)/dr \ dR F"'{R) f dR 
r ^ 2B{r) J ~d^ ^ F"{R) \Jh 



(47) 



0, 



dajr) 1 2 F"{R) dR\ F"{R)d^R F'" {R) 
dr \r ^ F'{R) J ^ F'{R) ^ F'{R) 



dR 
dr 



= 0. 



(48) 



Once F{R) is given, together with Ea. (|^5|) . the above equations form a system of three differential 
equations in the quantities a(r), B(r) and R{r). 

As already said, the static solutions describe a black hole if there exists a real positive solution 
vh of B(rH) = 0,B'{rH) 0. If this happens th is the radius of the event horizon. The Killing 
surface gravity reads 



d{e^o.(r)B{r)) dB{r) 



dr 



dr 



iga(''ff) 



dB{r) 



dr 



(49) 



Non trivial examples of such F{R) gravity black hole solutions have been discussed in [T3l and we 
shall deal with them in the next Sections. In the following, we shall show that the First Law of 
black hole thermodynamics holds, making use of the equations of motion and of the expressions 
for the Hawking temperature and entropy obtained by independent methods. First, the tunneling 
method gives for the Hawking temperature 



Tk = 



e"('«) dBirn) 



An 



dr 



(50) 



The entropy associated to these black holes solutions can be calculated by the Wald method 
(see Appendix II). One has 

Sw^^F'iRn). (51) 

For simplicity we will consider only spherical horizons, thus the area is Ah = 'inr'j^ and the volume 
Vh = I "''"If • By evaluating the equation of motion (|47p on the event horizon, and multiplying 
both sides of equation by drn, we have 

T,dS^ = e"(-) (g - ''"^^" rj,^ dr, . (52) 

Thus, we have derived for a generic F{R) gravitational model the First Law of black hole thermo- 
dynamics as soon as the identification 

can be made. Within these F{R) modified gravity theories, this is one of the main result of this 
paper. In the next Sections, by making use of several exact solutions, we will provide a support 
for this identification. 

Our proposal, expressed by Eq. (|53L should be compared with a similar proposal contained 
in Ref. [33]. In Ref. [3S] an attempt to define a local Misner-Sharp mass has been presented. 
There, however, the proposed formula is not really satisfactory, because the quasi-local form is 
only present in some particular cases, one of which will be discussed in the next Section. 
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3.1 The constant curvature case 

As a simple but important example, let us consider the class of static solutions with constant 
curvature Rq . In this case one has the solution with a = (in I12p , and the the first equation of 
motion (Ii7|l reduces to 

±(r-rB{r) + ^] ^0, (54) 



dr \ ' ' 3 

where 



A.^^^^. (55) 



Thus, one arrives at Schwarzschild-de Sitter solution 

5(0= (1-7-Aoy) , (56) 
and i?o = 4Ao. Here C is a constant of integration. The horizon is located a.t r = vh, where 

l = ^ + Ao^. (57) 

TH 3 



Making use of Eq. (|53p one has 



and by Eg. (1571) one finally gets 



E/t = ^ ( Fi(r„ - ^) I , (58) 



which is our identification of mass-energy expression for this class of black hole, in agreement with 
Ref.m. 



4 The Clifton-Barrow solution 

Let us apply the same procedure for the highly non-trivial Clifton-Barrow solutionfT^, for which 
a is not a constant. The starting point is the following F(i?)-modified gravity model: 

F{R) = R^+^G^ . (60) 

For dimensional reasons we have also included the Newton constant , 5 being a numerical 
parameter. When 6 = the Hilbert-Einstein action of GR is recovered. Note that in this case the 
modification with respect GR is not additive. 

Looking for static, spherically symmetric metric of the type ([T^ . we find the Clifton-Barrow 
solution of Eqs. (|in])-(|15]), that it 



i(r) ^ ±_ 



r \ *(i+2*)/(i-*) f{l~2S + 452)(i - 26- 26^) 



2 



(61) 



(1 — (5)2 / C \ 

^^'^^ ^ {1-26 + 452)(i -26- 252) - ^(1-25+45^)7(1-5) J ' (62) 

and 

cs 66(1 + 6) 

^=^' ^ (252 -f 25-1) (63) 

Above, ro > is an arbitrary constant while C > is the integration constant of the model. We 
assume 5^1. 
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The horizon radius, defined by B{rH) = and drB{rH) 7^ reads 



TH 



and since C > 0, the Clifton-Barrow metric is a black hole solution. 
According to Equation the Killing-horizon surface gravity reads 



" 2Y (1-2(5-2,52) ^5(1+25)7(1-5) ' 

which can be used to find the Killing- Hawking temperature Tk = kk/Stt. 

With regard to the black hole entropy associated with the event horizon of the Clifton-Barrow 
solution, from the Wald formula in Equation ([5T|) we find [46]: 



(2(52 + 25 - 



(66) 



In order to have the positive sign of entropy, we must require S > {VS — l)/2 or — 1 < (5 < 0. 
The solutions with < S < {VS — l)/2 or (5 < —1 are unphysical, whereas for 5 = we find the 
result of General Relativity. On the other hand, only the solutions of — 1 < (5 < give a real 
value for the Killing surface gravity kh- If <5 > (\/3 — l)/2 the Hawking Temperature becomes 
imaginary. 

Making use of Eqs. ([5^ one has 

dEK=Asr''f-'^/^'-'UrH, As = ^e''" G'-\cs)' [2(1 + 5) ~ cs 5] . (67) 
As a result, the energy turns out to be 

Finally, from Eq. (|64p one gets again that the energy is proportional to the constant of integration 
of the BH solution since 

= 5(1+25)7(1-5) ^ ' (69) 
where we have introduced the dimensionless constant depending on 5 

^ / 2^-W((5-l)2((5+ 1)^+1 1 \ 

' Wl-2(5-2^2^1_2,5 + 452 (2(52 + 2(5 - 1)V • ^ > 

We conclude this Section with some remarks. In the above expression, the range of parameter 
(5 has to be restricted to the ranges already discussed in order to have a positive temperature 
and entropy. As a check, it is easy to show that in the limit S ^ 0, one gets the GR value 
C — 2EG. Furthermore, the Killing energy Ek and the Killing temperature depend on the 
dimensional constant tq, and we may take it proportional to Planck length VG- 

5 1/R Model 

As a further non trivial example, let us consider the following i^(_R)-model: 
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where h and 7 are positive, dimensional, arbitrary constants (we may choose, for example 7 — 
1/G^). In Ref.^lSj it has been shown that this model admits a static, spherically symmetric 
solution of the type 

/ \ 1/2 

.-<••) =(i) , ,72) 

and R — 6/{hr). Here vq is an arbitrary constant, which is present for dimensional reasons, and 
C is the integration constant. Let us consider the solution of B{rH) = 0, namely 

If we assume C > 0, it is easy to show that there exists always a simple zero rn > 0, which defines 
the event horizon, and so the above solution represents a black hole. With regard to the related 
entropy, Ea. (j5ip gives 

the entropy being positive, since 7 > 0. The Killing temperature associates with the horizon reads 



27r \6h ^ 2r% J \ro 

By computing the Killing energy from (j53p we have 

/i7 /iV^^,/ 9/2 6 7/2\ _ /i7 / 1 / 9/2 6 7/2 



Thus, making use of Ea. ([71| one arrives at 

E,,^(LV"c. (78) 



63G \ro 

Also in this case we can identify the integration constant of the model as a quantity proportional 
to the black hole Killing energy. 



6 The Deser-Sarioglu-Tekin topological black hole solutions 

In this Section, first we generalize the modified gravity black hole solution of Deser et al. [TT] , and 
then we shall show that also for these solutions the First Law of black hole thermodynamics is 
valid and the constant of integration is proportional to the Killing energy. 

For the sake of simplicity we shall restrict ourselves to the four-dimensional case, but, since 
we are interested in black hole with generalized topological horizon, we have to include a non 
vanishing cosmological constant (see for example the GR case [TTl H51 The D-dimensional 

case as well as the inclusion of Electromagnetism presents no difficulties. 

To begin with, we write down the action of the model 

/ = / d'^x (r-2A + VSa Vf) , (79) 
lOTT Jm ^ ' 

where cr is a real dimensionless parameter and F = C^i/piC'"^'"' is the square of the Weyl tensor. 
For (7 = the Weyl contribution turns off and GR result is recovered. This model is a very 
interesting additive modification of GR with cosmological constant. 
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For more generality we look for static, (pseudo)-spherically symmetric solutions with various 
topology and so we write the metric in the form 



ds' = -a\r)B{r)dt' + + r' + P'd<l>'^ - (80) 

where the horizon manifold will be a sphere 52, a torus T2 or a compact hyperbolic manifold 
according to whether k = 1, 0, —1. 

A direct computation shows that the noteworthy properties of the Weyl scalar F discussed in 
Ref.fll] for fc = 1, are still valid for fc = 0, —1. Thus the unknown functions a{r) and B{r) can be 
obtained by imposing the stationary condition 51 = 0, where, / is the original action evaluated on 
the metric (j80p (up to integration by parts and on the "topological" variable p, cf)) . It reads 

i^j dr {{I- (7)[ra'{r)B{r) +ka{r)]+iaa{r)B{r) - Kr'^a{r)] , (81) 

from which it follows 

(l-(T)m'(r)+3CTa(r) =0, (82) 



f 1 - 4ct") 

rB'(r) + ^ "-lB{r) = k-K . (83) 

1 — cr 1 — (T 



Here we are assuming a ^ 1, ^ |. The general solutions are 



r 



air) = ( - ) , (84) 



^(^)-^(r^-^^"^-^3(i^' 

C and To being integration constants. 

One can see that black hole solutions exists only for negative cosmological constant, but in 
the case fc = 1, already discussed in |46| . where A can assume any arbitrary value. As usual, the 
horizon is given by the positive root rn of B(r) = with B' [rn) 7^ 0. The algebraic equation can 
be easily solved for the integration constant C and gives 



1 — (T . r 



2 



l-4cr 



C=[k A , ^ J rJr- . (86) 

V 1-4(7 3(l-2cr)y ^ ^ ' 

The equation (|83p evaluated on the horizon leads to 

THB'{TH)^k~K^. (87) 

1 — (T 

Thus, since the Killing-Hawking temperature Tr- = aHB'{rH)/4:n, taking into account Eq. (|84p 
we get 



ATTrnTK = fc - A 



{I- a)) yro 

On the other hand, a direct computation along the line discussed in [35] for the case fc = 1 leads 
to the BH entropy 

Sw^^(l + <y) , (89) 
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where, in order to deal with a positive entropy, we have to restrict to the interval a G (—1, !)■ 
Above, Ah ~ ^tr'jj, in which Vi = An (the sphere), Vq = |9t|, with t the Teichmueller parameter 
for the torus, and finally V-i — 'ing, g > 2, for the compact hyperbolic manifold with genus g 

As a result we have 
Furthermore, Eq. gives 

dC = {k- K^^^^^-^r^drn ■ (91) 
As a consequence the first Law holds and 

(», 

In this class of modified gravitational models the energy of black hole is particularly simple, since 
the modification is describe by the dimensionless parameter a. 



7 Topological Conformal Weyl Gravity 



In this Section, first we revisit the higher gravity black hole solution of Riegert and others |50l[5T] 
and its topological version [52J . 

To begin with, we write down the action of the model in the form 



/= / d'^x^/^["i[R~2K) + ZuF] 

M 



(93) 



where 7 is an arbitrary parameter, which may be proportional to the square of Plank mass, uj 
is a dimensionless parameter and F — C^^psC^'^^^ is the square of the Weyl tensor. The pure 
conformal invariant model 7 = is very interesting and its phenomenology has been investigated 
in Ref.p,. 

As in previous Section, also here we shall consider various topology and this means that the 
metric will have the form ([5D|) . and the arbitrary functions a{r),B{r) will be obtained from the 
reduced action 



I = dr 



7 [rB{r)a'{r) + ka{r) - 2A r^a(r)) + w 



A^{r) 
r^a{r) 



(94) 



where we have put 
Mr) 



r^a{r)B"{r) + ir'^a {r)B' {r) - 2ra{r)B'{r) + 2r^a"{r)B{r) 
-2ra'{r)B{r) + 2a{r)B{r) - 2ka{r) . 



(95) 



As a result, we are dealing with a higher order Lagrangian system, the Lagrangian depending on 
the first and second derivative of the unknown functions a(r) and B{r). 
The equations of motion read 



dr"^ 



AB{r) 
a{r) 
A 

I 



A 



dr \ra{r) 



[3rB'{r) - 2B{r)' 



[r'^B"{r) - 2rB'{r) + 2B{r) - 2k\ 



A^ 7 

h — 

r'^a?-{r) lo 



[k - B{r) - rB'{r) - 2Kr'^] = , 



(96) 
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(fA(r) d ( Air) , 

^ ' I ^ -[ira'ir) - 2a{r)] 



dr^ dr \ra{r) 

+ 4t\ [2r^«"(0 - 2ra'(r) + 2a{r)] + = , (97) 

r'^a{r) Zuj 

For simplicity let us look for exact solutions with a(r) = 1. With this Ansatz Eq. ((97)1 can be 
integrated and one obtains 

B{r) = - + Co + cir + csr^ , (98) 
r 

hi and Ck being integration constants. In order to satisfy Eq. (j96p we have to distinguish the two 
cases 7 (a modified Einstein gravity) and 7 = 0, (pure conformal gravity), since they provide 
completely different solutions. 

In the case 7 7^ Eq. (|96p is satisfied only if 

co = fc, ci = C2 = C3 = -iA, (99) 

while hi remains a free parameter. We see that this is a topological Schwarzschild-de Sitter(AdS) 
black hole like solution, since 

B(r)=fc---iAr2, (100) 
r 6 

where here hi has been replaced by C . It has to be noted that this is the solution which one would 
have obtained from the Hilbert-Einstein action with cosmological constant, that is with u = 0. 

As we already said, if 7 = the solution is completely different and in fact, in such a case 
Eq. is satisfied only if 

„2 _ t,2 

Now the solution depends on the three arbitrary parameters cq, C2 and hi. By a redefinition of 
them by Co — > fc + 3co, C2 — >■ A, fei — > — C, we write it in the form 

B{r) =k + icQ-^ (2fc + 3co) r + Ar^ - - , (102) 
G r 

in agreement with the topological black hole solution already found by Klemm in |52j. 

The event horizon exists as soon as there is positive solution rn of B{r) — 0. For example, if 
C > and A = > 0, it is easy to show that there exists always a positive root independently 
on the values of co and of L, while, in the opposite case A < 0, a positive root of B{r) = exists 
only if Co > and the value of |A| is sufficiently small. The special A = case will be discussed at 
the end of this Section. 

With regard to the computation of Entropy, assuming that there exists an event horizon 
B{rH) = 0, with rn > and B'{rH) ^ 0, for the pure Weyl gravity case the Wald method gives 

Sw = 2ujVk ( Co I = 2ujVk (x — co) dSw = 2ujVk dx , (103) 

Vh J 

where for convenience we have introduced the variable x = C /rn- Here Ah = VkT% {k = 1,0,-1), 
with Vi = 47r, for the sphere, Vb = |3t|, r being the Teichmueller parameter for the torus, and 
V-i — 4:Trg, g > 2, for the compact hyperbolic manifold with genus g [47| . The integration constant 
C in Eq. p03p can be seen as a function of rn obtained by solving the equation Birn) — 0, which, 
as it follows from (|102p . it is a second-order algebraic equation in C. Of course, in order to have 
a positive entropy we have to choose co < C /rn ~ x and moreover C has to be positive being 
proportional to the energy. 
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Now we restrict ourselves to the A = l/L^ > case. In this way, by solving the equation 
Birn) — with respect to C we get 



2r* \^ 
2x= — --| + fc + 3co + Vw^, W^= -| + fc + 3co - 4co(2fc + 3co) > , (104) 
TH ) 

and from the latter equation it follows 

da; -^i- 1 + ^ = drn , dC = rudx + xdrH. (105) 

y VW J 

On the other hand the Hawking temperature can be written in the convenient form 

B'{r) 



Tk = 



H 1 C^rl 



47r 47r th \ L'^ 

and using Eqs. p03l) and (I105|) we obtain 



W , (106) 



ttL 



^{rndx + X drn) ■ (107) 



We finally see that the First Law of black hole thermodynamics reads 



TKdSw = '^dC . (108) 

As a result, we may again identify the energy as 

Ek = '^C. (109) 

We conclude this Section with some remarks. The pure Weyl conformal gravity does not contain 
dimensional parameters. Thus, one could think that there exists a trivial entropy and a vanishing 
energy, but, as we have shown above, the solution gives rise to a length scale L related to the 
integration constant A. In such a case the First Law of black hole thermodynamics holds and the 
energy of black hole solution is proportional to the other dimensional constant of integration C. 

The situation is different when A = 0, since in such a case the scale does not emerge and for 
the horizon one gets 



rH _ + 3co + ^(cq + fc)(fc - 3co) 

C ~ 2co(2/c + 3co) ■ ^ ' 

The latter equation gives a positive rn for k ^ and a suitable value for cq. In any case we see 
that X = C/rn is a pure number and so da: = and the entropy is trivially constant. The First 
Law of black hole thermodynamics is trivially valid with a vanishing energy. This is the particular 
case discussed in [52]. 



8 Conclusions 

In this paper the issue of defining the energy associated with a static, spherically symmetric black 
hole solution in higher order modified gravitational models has been tackled. We have proposed 
to identify the black hole energy as a quantity proportional to the constant of integration, which 
appears in the explicit black hole solution. The identification is substantiated by the fact that 
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in all explicit and known examples, we have been able to show that the First Law of black 
hole thermodynamics (Clausius relation) holds true as a consequence of equations of motion, 
and evaluating in an independent way the related entropy via Wald method and the Killing- 
Hawking temperature via quantum mechanics techniques in curved space time. In the case of 
F{R) modified gravity some non trivial exact black hole solutions have been considered. In the 
case of modified gravity in which the quadratic Weyl scalar is additively present, first we have 
found the corresponding new topological black hole solution, and then we have verified for it 
our proposal. Finally our proposal has been shown to work also in another non trivial higher 
order gravity theory, namely the topological conformal Weyl gravity. It is easy to show that the 
proposal is also working for constant curvature black holes solutions in the usual Einstein gravity 
with cosmological constant modified by generic curvature-squared terms (see, for example [54]). 

On general grounds, we may say that our explicit results are in agreement with the general 
result obtained recently in Ref. [28 , and together other results appeared in literature seem to in- 
dicate that the thermodynamic origin of a generalized modified gravity, when horizons are present, 
has a broad validity. 



9 Appendix I: the Tunneling method 

In this Appendix, for the sake of completeness, we present a short review of the tunneling method 
in its Hamilton- Jacobi variant [25]. The method is based on the computation of the classical action 
/ along a trajectory starting slightly behind the trapping horizon but ending in the bulk, and the 
associated WKB approximation (c = 1) 

Amplitude cx 6*"^ . (Ill) 

The related semi-classical emission rate reads 

r (x lAmplitudep cx e"^^^^ . (112) 

The imaginary part of the classical action is due to deformation of the integration path according 
to the Feynman prescription, in order to avoid the divergence present on the horizon. As a result, 
one asymptotically gets a Boltzmann factor, in which an energy w appears, i.e. 

Fcxe"^", (113) 

and the Hawking temperature is T = -g. 

To evaluate the action /, let us start with a generic static, spherically symmetric solution in 
D-dimension, written in Eddington-Finkelstein gauge, which, as it is well known, is regular gauge 
on the horizon 

ds^ = -B{r)e^°'^''Uv^ + 2e"^'''>drdv + r^dnl,_2 = ^ij{x')dx^dx^ + dn\,_^ (114) 

here = (u, r), where v is the advanced time. Since we are dealing with static, spherically sym- 
metric solution space-times, one may restrict to radial trajectories, and only the two-dimensional 
normal metric is relevant, and the Hamilton- Jacobi equation for a (massless) particle is 



-f'^djdjl = +2e"(-'-'>dJdrI + e2"('')B(r)(a^/)2 = . (115) 

Thus 

in which uj = —dyl is the Killing energy of the emitted particle. In the near horizon approximation, 
B{r) ~ B'^{r ~ rn)- As a consequence, making use of Feynman prescription for the simple pole 
in r — r/f . one has 

drdrl - / dr———^ , (117) 
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where the range of integration over r contains the location of the horizon rn ■ Thus 



^/=J^, (118) 



and the Hawking-KiUing temperature is 



T. = l^. (119) 

If one had introduced the Kodama energy ujh — e~°'"uj, one would have obtained the Hayward 
temperature Tjj — ■ 



10 Appendix II: the Wald Entropy 

In this Appendix, we recall the basic formula for computing black hole entropy for a generalized 
covariant theory of gravity. Following [271 [551 [55] . the explicit calculation of the black hole entropy 
Sw is provided by the formula 

Sw ^ ~'2Tr [ ( Y^—] eape^sdS , (120) 



S X^RaP'jS J }j 

where ^ — ^(Ras-ys) is the Lagrangian density of any general theory of gravity, = —epa, is 
the binormal vector to the (bifurcate) horizon. It is normalized so that Cape"^ = —2. 
For the metric ([T^ . the binormal turns out to be 

eo/5 = e"M(<5°<5i-5i^^). (121) 

The induced area form, on the bifurcate surface r = rn, t ^constant, is represented by dS. 
Finally, the subscript (iJ) indicates that the partial derivative is evaluated on the horizon, and 
the variation of the Lagrangian density with respect to Ra^-yS is performed as if Rap-yS and the 
metric gaii are independent. For example, 



5R 



- ^ {g'^^g"^' - g'^-g^^) ■ (122) 



^Rlivap 2 

As a result, for the modified gravity model of the F{R) class, one obtains 
This is the formula used in Section III and IV. 



(123) 
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